The cross-correlation search for a hot spot of gravitational waves 
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The cross-correlation search has been previously applied to map the gravitational wave (GW) 
stochastic background in the sky and also to target GW from rotating neutron stars/pulsars. Here 
we investigate how the cross-correlation method can be used to target a small region in the sky 
spanning at most a few pixels, where a pixel in the sky is determined by the diffraction limit which 
depends on the (i) baseline joining a pair of detectors and (ii) detector bandwidth. Ifere as one of 
the promising targets, we consider the Virgo cluster - a "hot spot" spanning few pixels - which could 
contain, as estimates suggest ~ 10^^ neutron stars, of which a small fraction would continuously emit 
GW in the bandwidth of the detectors. For the detector baselines, we consider advanced detector 
pairs among LCGT, LIGO, Virgo, ET etc. Our results show that sufficient signal to noise can be 
accumulated with integration times of the order of a year. The results improve for the multibaseline 
search. This analysis could as well be applied to other likely hot spots in the sky and other possible 
pairs of detectors. 

PACS numbers: 95.85.Sz,04.80.Nn,07.05.Kf,95.55.Ym 



I. INTRODUCTION 

An enigmatic prediction of Einstein's general theory 
of relativity are gravitational waves (GW). With the ob- 
served decay in the orbit of the Hulse- Taylor binary pul- 
sar agreeing within a fraction of a percent with the theo- 
retically computed decay from Einstein's theory, the exis- 
tence of GW was firmly established. Currently there is a 
worldwide effort to detect GW with the operating inter- 
ferometric gravitational wave observatories, the LIGO, 
Virgo, GEO and TAMA [J. Now the advanced detec- 
tors being constructed include the upgraded LIGO and 
Virgo, the LGGT of Japan, LIGO-Australia and future 
possibilities such as Einstein Telescope (ET) [2]. 

Different types of GW sources have been predicted and 
may be directly observed by these advanced detectors in 
the near future (see 3] and references therein for recent 
reviews). In this paper we will address the problem of 
the targeted search of stochastic GW from a small region 
in the sky, typically of linear size of a few degrees (few 
pixels - a pixel determined by the diffraction limit) - a 
"hot spot" - where there is likely to be an abundance of 
independent, unresolved GW sources continuously pro- 
ducing a relatively large stochastic background. Such a 
scenario seems feasible for the Virgo cluster, which could 
contain about 10^^ neutron stars, the current estimate 
being 10* — 10® per galaxy. Out of these neutron stars 
a small fraction of them could be rotating sufficiently 
rapidly emitting GW in the advanced detector bandwidth 
of several 100 Hz to about 1 kHz. These could produce 



a reasonable signal-to-noise ratio (SNR) with an integra- 
tion time of the order of an year. Thus, the GW source 
consists of spinning asymmetric neutron stars whose am- 
plitudes and phases are randomly distributed. We will 
be thus dealing with a localized stochastic GW source. 
This is only one type of GW source, but there could be 
contributions from other sources such as supernovae with 
asymmetric core collapse, binary black hole mergers, low- 
mass X-ray binaries and hydrodynamical instabilities in 
neutron stars, or even GWs from astrophysical objects 
that we never knew existed. These will only in general 
(statistically) add to the SNR. The detectors we consider 
for this paper are advanced detectors such as the LIGO, 
Virgo, LGGT, ET etc. which are expected to have suffi- 
cient sensitivity for detecting a hot spot. 

The appropriate method for observing such a source 
is the cross-correlation method described in Q (hence- 
forth referred to as paper I), which is also generally 
known as the radiometric method. The idea is to cross- 
correlate data streams from two detectors with an ap- 
propriate time-delay, namely, the time-delay between ar- 
rival times of a GW wavefront from a specific direction 
Cl. This choice of time-delay allows the sampling of the 
same wavefront. As the detector baseline rotates with the 
earth, the time-delay between the data streams changes 
during the course of the day. The statistic targets a patch 
(pixel) in the sky around Cl its size being determined by 
the diffraction limit, namely, the inverse of the band- 
width divided by the light travel time along the baseline. 
This statistic in fact is a point estimate of the signal re- 
ceived from the given direction Q, and is most appropriate 
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for observing a hot spot and could be made optimal by 
'masking' the rest of the sky if the hot spot emits a strong 
signal. 

The GW strain amplitude for a rotating neutron star 
is proportional to the square of the frequency |5|] , 

/^-4.^«^|/^ (1) 

where a < 1 is the orientation factor, G is the Newton's 
gravitational constant, c the speed of light, e is the el- 
lipticity of the neutron star, I the moment of inertia, 
R the distance to the source and / the GW frequency. 
Since the cross-correlation statistic is quadratic in the 
strain amplitude, it scales as the fourth power of the fre- 
quency and therefore the main contribution to the SNR 
will tend to come from high frequency sources assuming 
that they are relatively abundant in the high frequency 
regime. Thus it is the population of millisecond neutron 
stars that we must primarily consider. We then estimate 
the millisecond neutron star population from the astro- 
physical information that is available and show that one 
can get an acceptable SNR, p ~ 1, for an integration of 
time of about an year. Using multiple baselines improves 
the SNR further. We find that among the current or near 
future baselines, the baseline of the two LIGOs and the 
baseline of LIGO Livingston and a LIGO like detector at 
AIGO site stand out - they give dominant contribution 
to the SNR. 

In section we give a brief description of the cross- 
correlation method and the statistic and then derive an 
expression for the optimal SNR. In section Hill we state 
our results and discuss them in light of the astrophysical 
scenarios that are possible and the sensitivities of the 
future advanced detectors such as the ET. 



II. THE CROSS-CORRELATION STATISTIC 
FOR TARGETING A HOT SPOT 

We refer to paper I for the detailed arguments involved 
in defining the cross-correlation statistic. Here we only 
furnish the salient steps. Since here we are interested in 
observing a hot spot, we will restrict our discussion to a 
point source. The full statistic, which we denote by S, is 
a weighted sum of elementary pieces A^fc, k = 1,2, ...n 
defined over a time-segments — At/2 < t < tk + ^t/2 
which are labeled by k. The full observation time is T = 
n At. The At is so chosen that it is much larger than the 
possible time-delay between the detectors (which must 
be less than about 40 ms for ground-based detectors) 
and much less than the time required for the orientation 
of the detectors to change appreciably and also on the 
timescale in which the noise is stationary. Current values 
of At used in LSC data analysis vary from 32 to 192 
seconds. Let us consider a pair of detectors labeled by 
/ ~ 1,2, then the data in the I^^ detector is given by 
xi{t) = hi{t) + ni{t), the signal hj{t) is added to the 
noise ni{t) in the I^^ detector. For a point source in the 



direction fl, the ASk also becomes a function of Cl. It 
can be expressed easily in the Fourier domain, 

/oo 
dfIUtk;f)x2{tk;f)Q{tkJ,Cl), (2) 
-oo 

where the x*j(tk]f) are short term Fourier transforms 
(SFT) defined only over the interval At around tk, 
namely, 

xi{tk;f) / dt'xi{t')e-^^'f' . (3) 

Jtk-At/2 

The Q{tk, f,Cl) is a filter function chosen so that it op- 
timizes the filter output. It also depends on the power 
spectrum of the GW source and the power spectral den- 
sities of the noises in each of the detectors. As discussed 
in paper I, in the general case it is a far more complicated 
object - a functional - but for the case of a point source, 
it reduces to a function of the direction ft. Even then 
it remains a functional of the signal power spectral den- 
sity and the noise power spectral density (PSD). With a 
slight abuse of notation we still write it as a function of 
/• 

The ASk are random variables because of the noise 
and for different k we take them to be uncorrelated. The 
mean and the variance of ASk are denoted respectively 
by Hk = {ASk) and = {ASD - {ASk)^- It has been 
shown in paper I that the linear combination that yields 
the maximum SNR is: 

P=|E/^?>^|\ (5) 

where p is the SNR. The sum over k can be converted 
into an integral over t and henceforth in this article we 
drop the suffix k and replace tk by just t. This helps to 
avoid clutter without jeopardizing clarity. 

We now turn to the noise and signal PSDs in terms 
of which the SNR can be finally expressed. The signal 
cross-correlation in the two detectors in the limit of large 
time segment can be written as: 

{hUtj)h2itJ')) = s{f-r)H{f)^it,f,n), (6) 

where ^{t, f,Ct) is the so called directed overlap reduction 
function analogous to the one defined in '6] for the full 
sky, and given in the case of the point source by, 

7(i,/,J^) = r(i,J^) e2"*^"-^"(*)/", (7) 

r(fj,t) - F+i{t,n)F+2{fi,t) 

+ Fxi(jl,t)Fx2(rj,t), (8) 

and where the Ax(t) is the vector joining detector 1 to 
detector 2 and rotates with the Earth tracing out a cone. 
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The F+i,Fxi, I = 1,2, are the antenna pattern func- 
tions for the two detectors and for the two polarizations. 
As mentioned in paper I the directed overlap reduction 
function has a bandwidth of about 750 Hz as compared 
to the few tens of Hz for the overlap reduction function 
found by integrating over the full sky. This is the main 
advantage of this method in which the sensitive region 
of the detector bandwidth is sampled by the statistic. 
Further the quantity pH{f) is essentially the flux per 
unit frequency per unit solid angle. For the noise, we 
take the noise in the two detectors to be uncorrelated, 
{ni{t)n2{t')) — 0, and the one-sided noise PSD in each 
detector / is given through the defining equation. 



{n}{t;f)nj{t;n) = ^S{f 



nPi{t;\f\)- (9) 



We also assume {hi{f)nj{f')) — 0, I,J — 1,2, that is 
the signal and noise are uncorrelated. We are now ready 
to write down the optimal filter. In paper I it has been 
shown that the optimal filter for a given time segment 
labeled by t and for a point source in the direction is 
given by. 



QitJ,Cl) = A(t) 



H{f)r{tj,n) 

P,{t;\f\)P2{t;\f\)-^ 



(10) 



where X{t) is a normalization constant, which in any case 
cancels out in the SNR. The SNR p is given in terms of 
fi{t) and a(t) which are the mean and standard deviation 
respectively of AS{t). To keep the expressions simple we 
assume that the noise in the detectors is stationary. This 
certainly will not be the case, but since we are only inter- 
ested in order of magnitude results, the assumption is not 
unjustified. Then Pj becomes a function of / only. Also 
we consider a band-width /i < / < /2 for evaluating the 
SNR; the lower limit /i is determined by the seismic cut- 
off, while the upper limit /2 is decided by the GW sources 
above which we do not expect significant contribution to 
the SNR. Given this, the relevant quantities can be best 
expressed in terms of the following two averages: 



>BW 



2 /"^^ 

A/ 4 Pi(/)F2(/)' 



(11) 



(r2)lday(0) - 



1 



Tl day Jo 



T^{n,t)dt, (12) 



where A/ = ji — j\- The first is the noise weighted av- 
erage of the signal H'^{f), the suffix BW denotes band- 
width, while the second is the time average of the squared 
directed overlap reduction function taken over one side- 
real day. It is a function of sky position of the source. 
But since the azimuth is averaged over 27r, it is just a 
function of the declination of the source. Then in terms 
of these averages we have. 



^I{t) - A (AtA/) {h^)bw r^{n,t) 
1 



ait) 



-A (AtAfY^ 



(13) 
(14) 





LIGO-H 


Virgo 


LCGT 


AIGO 


LIGO-L 


0.387 


0.288 


0.224 


0.452 


LIGO-H 




0.214 


0.215 


0.312 


Virgo 






0.276 


0.286 


LCGT 








0.256 



TABLE I: The values of the square root of the one day (side- 
real) average of are given for the Virgo cluster whose dec- 
lination is ~ +12.7° (the RA is irrelevant since we take a 
one day average). LIGO-L stands for LIGO-Livingston and 
LIGO-H for LIGO-Hanford. 



Then using the continuous limit of Eq. ([5]) , we may write 
the SNR p in terms of the averages as follows: 



1 

At 



dt 



pHt) 



2 (rA/)5 {H- 



>l/2 
/BW 



,1/2 
' 1 day ■ 



(15) 



We now use this expression to compute the SNR for the 
continuous wave sources from the Virgo cluster. We ob- 
serve that the SNR scales as ^/T. 

To fix ideas we can look at a simplified situation of 
identical detectors with white noise Pi{f) = Pq in the 
frequency range /i < / < /2 and Pj = oo otherwise. 
Similarly we may consider flat signal spectrum H[f) = 
Hq, then the SNR simplifies to: 



p = 2 [TA/]^ 



Ho 



,1/2 
' 1 day ■ 



(16) 



The values of (F^ 



1/2 

)i day various combinations of de- 
tector baselines are given in Table I for the Virgo cluster 
which has a declination ^ +12.7°. 



III. RESULTS AND DISCUSSION 

A. Pulsar population and distribution 

We consider gravitational waves from rotating neutron 
stars in the Virgo cluster. One important parameter 
of this source is the population of such neutron stars. 
The number of Galactic neutron stars is estimated to 
be 10^ — 10^ since the birth rate is about 10~^/yr and 
the age of the Galactic disk is about 10^°yr. What is 
more important in our case is the number of Galactic 
neutron stars whose rotation period is of the order of 
milliseconds. From the survey of radio pulsars in our 
Galactic disk, the population of millisecond pulsars is es- 
timated to be at least 40000 which implies a birth 
rate of 2.9 x 10~^/yr. This is consistent with other stud- 
ies of the millisecond pulsar population by Ferrario and 
Wickramasinghe (3.2 x 10~^/yr) [l^ and by Story et al. 
(4 — 5 X 10^^/yr) JT,]. From the recent observation of 
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FIG. 1: The distribution of observed radio pulsars. Tiie hor- 
izontal axis is logjQ(/r) where fr is the rotational frequency 
of pulsars. The histogram is the observed number. The solid 
line is the two component Gaussian model of the distribution. 



cluster. We assume that the total number of neutron 
stars in our Galaxy is 10® for < 50Hz, and 40000 for 
/,■ > 50Hz. Since there are approximately 10^ galaxies in 
the Virgo cluster, total number of neutron stars in Virgo 
cluster is Mow ~ 10" for fr < 50Hz, iVhigh ~ 4 x 10^ for 
/,■ > 50Hz. The distribution of neutron stars including 
millisecond pulsars in Virgo cluster thus becomes 

Nif)df = 

(iVMghPl(l0g/.) + ^l°w^'2(log/.))-^r^. 

(19) 

Since the length of data of one time-segment, At is at 
most 10'^ seconds, the frequency resolution is larger than 
lO"'^ Hz. The frequency bandwidth can be taken as 10'^ 
Hz. Thus the number of frequency bins is 10^. Since the 
number of pulsars with / > 100 Hz is 10^, the number 
of pulsars in each frequency bin is about 10. In the low 
frequency regime this number is much larger. Thus, it is 
not possible to resolve the signal from each pulsar, which 
confirms the stochastic nature of the Virgo cluster hot 
spot. 



gamma rays with the Fermi satellite the popula- 

tion of millisecond pulsars in our Galactic globular clus- 
ter is estimated to be 2600-4700 which is one order lower 
than millisecond pulsars in the Galactic disk. Although 
there might be significant population of millisecond pul- 
sars which do not emit radio waves, X and gamma rays 
now, since the life times of millisecond pulsars are be- 
lieved to be long 10-^°yr) [l^, we do not expect a large 
population of such millisecond pulsars to exist. Thus we 
adopt 40000 as a typical number of neutron stars per 
galaxy whose rotation period is of the order of millisec- 
onds. 

A catalog of radio pulsars is given in the ATNF pulsar 
database [1J| . The distribution of observed radio pulsars 
is given in FigHJ We find that the distribution naturally 
falls into two regions separated by 50Hz. In each region, 
the distribution is approximately Gaussian as seen from 
the figure. This means that the distributions in each 
region may be approximated as log-normal distributions 
given by: 



Pi(log/,)d(log/,) 



P2(l0g/.)d(l0g/,) 



1 _ (log/r-logui)'' 

e d(log/.). 



27r(Ti 

(for fr > 50Hz) , (17) 

1 _ (log /^-IogM2)^ 

-e ^1 d(log/,). 



(for fr < 50Hz) , (18) 



where m = 219Hz, cti = 0.238, /i2 = 1.71Hz and 
0-2 = 0.420, and fr — //2 (/ is the gravitational wave 
frequency). Pi and P2 are normalized to unity when in- 
tegrated from = to infinity. We assume a similar bi- 
modal form of distribution of neutron stars in the Virgo 



B. Signal-to- noise ratio 

The spectral density of gravitational radiation from 
neutron stars in Virgo cluster, H{f) is given as 

Hif) 

= {h')Nif) 

7.05 X 10-34 f ^) 
x{a')fNif), 



1.1 X lO'^^gcm^ 



(20) 



where (a^) represents the average with respect to the 
inclination angle and the polarization angle. Assuming 
uniform distribution of the sources over the angles, we 
have (a^) — 0.4. We also used the distance R — 16.5Mpc. 

In order to obtain a rough idea of how large H{f) 
is compared with the noise power spectral density, it is 
convenient to define an effective source power, Hcsif) 

by, 



8Tobs(r^)l day/-ff^(/), 



(21) 



where Tobs is the observation time. 

Then the signal-to-noise ratio is given by. 



/2 



/l 



/ Pl{f)P2{f) 



1/2 



(22) 



The noise power spectral density of various advanced de- 
tectors including Einstein Telescope, as well as Hcs{f) 
are plotted in Figl2] Here, we assume Tobs — lyr and 
0.2. In this plot, Hce{f) is plotted for 



/I day 
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e = 10~^, 10~^, 10""^. Although in these plots, we in- 
clude the contribution from low frequency neutron stars 
with fr < 50Hz, the contribution of these to the SNR is 
very small, only about a few percent. 




FIG. 2: One-sided noise power spectral density of LCGT, ad- 
vanced LIGO, advanced Virgo, and Einstein Telescope (ET- 
B). LCGT noise curve is "variable RSE in broadband mode" 
(VRSE(B)) Advanced LIGO noise curve is "Zero Det, 

High Power" taken from [l^ . Advanced Virgo noise curve 
is take from the Virgo websi te IitI] . The Einstein Telescope 
noise curve is called "ET-B" [la ]. The effective source power 
^cS^if) is S'lso plotted. In this plot, we assume Tobs = lyr 
and {FX/' = 0.2. 



We now consider the quantity {H'^)^^Af^/'^. We find 
that, 



cx 



VIO-V ll.lxlO' 



lO^SgcmV V4xl0 



(23) 



where iVmsp is a numbers of millisecond pulsars. The 
signal-to-noise ratio with 1 year observation is written as 



/lobs / (-'^^)l^day\ / £ 



P = Pl.vrl 



lyr \ 0.2 / VlO 



(lO-s) 



1.1 X 1045gcmV V4xl07 



msp 



(24) 



we take p > 3, there is more than 99.7% chance that the 
noise is not masquerading as the signal. We then have 
the following result: 




(25) 



The values of piyr and T^^^^ are given in Tables II and 
III for (r2)J/^^y = 0.2. These tables along with Eq.^ 
and Eg. ipS)) can be used to obtain the p and the Tobs for 

any other value of (r^)^^^^^^ or equivalently for any other 
sky location, not just the Virgo cluster. From the tables, 
we find that in the case of e = 10~^ we can achieve p — 3 
in about 3 months with advanced LIGO noise PSD. For 
advanced Virgo and LCGT, it takes about 1.5 year to 
achieve it. For Einstein Telescope, it will be quite easy 
to observe it. However, the results strongly depend on 
the value of e. If e = 10^^, it will become difficult to 
observe the Virgo cluster hot spot with advanced LIGO, 
advanced Virgo and LCGT. Only Einstein Telescope will 
be able to detect it. Note that these are order of mag- 
nitude results where we have assumed a typical value of 

(r^>lly-o.2. 

In Table IV and V, the values of piyr and T^^^'^ ob- 
tained for various detector combinations of advanced de- 
tectors, are given for the specific source location of the 
Virgo cluster. Since (r^)y^^y is roughly factor of 2 larger 
than 0.2 for the LIGOs and AIGO network, the values of 
piyr and T^j^'^ are improved significantly for these base- 
lines. For example, the piyr = 3 is achieved in 26 days 
by LIGO-L and LIGO-H, and in 19 days by LIGO-L and 
AIGO. 



Plyv 


LIGO 


Virgo 


LCGT 


ET-B 


LIGO 


5.8 


2.5 


3.0 


52 


Virgo 




1.5 


1.4 


24 


LCGT 






1.6 


27 


ET-B 








4.7 X 10^ 



TABLE II: The signal-to-noise ratio piyr which can be ob- 
tained with 1 year observation time for each combination 

of the detectors' noise PSD assuming (r^)^^ = 0.2 and 



£ = 10" 



We now compute the observation time required to achieve 
p — 3, which we denote by T^j^'^. We choose this value of 
SNR because the noise in the statistic S{d), as argued 
in paper I, is distributed as a Gaussian with mean ps 
and standard deviation erg. This is the consequence of 
the generalized central limit theorem. The p is ps/'^S- 
When no signal is present, we have, ps = ^- Thus when 



The results improve if we employ several baselines of a 
network of detectors. A full treatment of multi-baseline 
gravitational wave radiometry has been given in jigj . The 
results of this paper can be easily applied to the case of 
the hotspot where the source consists of a single pixel or 
at most a few pixels. In this case the beam matrix for 
a single baseline essentially consists of a single diagonal 
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rrip—'i 

-'-ohs 


LIGO 


Virgo 


LCGT 


ET-B 


LIGO 


0.26 yr 


1.5 yr 


1.0 yr 


1.2 day 


Virgo 




4.2 yr 


4.7 yr 


5.7 day 


LCGT 






3.6 yr 


4.5 day 


ET-B 








1.3 X lO^sec 



TABLE IIL Observation time T^^^^ required to achieve p = 
3 for each combination of the noise PSD and assuming 



day 



: 0.2 and e = 10" 



Plyr 


LIGO-H 


Virgo 


LCGT 


AIGO 


LIGO-L 


11.3 


3.54 


3.36 


13.2 


LIGO-H 




2.63 


3.21 


9.12 


Virgo 






1.90 


3.51 


LCGT 








3.84 



TABLE IV: The signal-to-noise ratio piyr which can be ob- 
tained witli 1 year observation time for each combination of 

the detectors' noise PSD and (P^)? day in Table |l] e = 10"^ 
is assumed. Noise PSD of AIGO is assumed to be the same 
as that of LIGO noise PSD. 



term for a single pixel or in case of few pixels, a small 
block diagonal matrix having dominant diagonal terms. 
The p which we have defined above is then just the SNR 
obtained for the log likelihood statistic A defined in that 
paper. We also deduce from further results of that paper 
on sensitivity that approximately in our case, 



2 _ \ " 2 

^network — ' 
I 



(26) 



where Pnctwork is the SNR for the network and the index 
T runs over all the baselines of the network. Similarly, 
it is easy from the foregoing to deduce that the obser- 
vation times to reach an SNR of 3, namely T^^^^ , add 
harmonically; more specifically we have. 



1 



rpp=3 

^ohs 



^ ^ rpPX=3 

I ^obs 



(27) 



where now the T^^^ denotes the time of observation re- 



Tolf [day] 


LIGO-H 


Virgo 


LCGT 


AIGO 


LIGO-L 


25.8 


262 


291 


18.9 


LIGO-H 




474 


319 


39.5 


Virgo 






907 


266 


LCGT 








223 



TABLE V: Observation time T^^^^ required to achieve p = 3 

for each combination of the noise PSD and (P'^)? day Table 
[T] e = 10~^ is assumed. Noise PSD of AIGO is assumed to 
be the same as that of LIGO noise PSD. 



Detector combination 


plyr 


r^bf [day] 


L-H-V 


12.1 


22.4 


L-H-J 


12.2 


22.1 


L-H-A 


19.6 


8.55 


L-V-J 


5.24 


120 


L-V-A 


14.1 


16.5 


L-J-A 


14.1 


16.4 


H-V-J 


4.57 


157 


L-V-A 


10.1 


32.1 


L-J-A 


10.4 


30.3 


V-J-A 


5.54 


107 


L-H-V-J 


13.1 


19.1 


L-H-V-A 


20.4 


7.90 


L-H-J-A 


20.5 


7.81 


L-V-J-A 


15.1 


14.4 


H-V-J-A 


11.5 


25.0 


L-H-V- J-A 


21.4 


7.21 



TABLE VL The signal-to-noise ratio piyr which can be ob- 
tained with 1 year observation time and the observation time 
required to achieve p = 3 by more than 2 detectors. These 
are derived from Eqs. (f26)l and (p7)l and Tables IV and V. 
L: LIGO-Livingston, H: LIGO-Hanford, V: Virgo, J: LCGT 
in Japan, A: a detector with LIGO's noise PSD at the AIGO 
site in Australia. 



quired for the network and T^^~^ denotes the observation 
time required for the baseline I to reach the SNR of 3. 

We can now apply these results to various networks. 
The results are given in Table VI. We first consider the 
3 detectors, LIGO- Virgo (LHV) network. Just compar- 
ing tables IV and VI, the piyr goes up from 11.3 for two 
LIGOs to 12.1 for L-H-V network which is about 7 % 
increase. Note that one must here take into account 3 
baselines L-H, L-V and H-V. The Tj^^^ comes down from 
25.8 days for the two LIGOs to 22.4 day for the L-H-V 
network which is a decrease of 13 %. If one considers the 
two LIGOs along with the LCGT the improvement is al- 
most similar to Virgo case, that is, the observation time 
comes down to 22.1 day. The improvement of adding 
other baselines to the L-H baseline is marginal because 
the LH contribution is dominant. Note however that an 
interesting improvement is obtained if one considers a 
detector at AIGO site assuming same noise PSD as the 
LIGOs. In such a L-H-A network, L-A contribution be- 

comes dominant because of largest (r'^)l in Table I, 

and Plyr goes up to 19.6 and T^^^ comes down to 8.5 
days. L-V-A and L-J-A networks are similar and gives 
second largest value of piyi- among 3 detector networks. 
They are better than L-H-V and L-H-J cases. 

In case of a 4 or 5 detector network, we can have further 
improvement, but the effect is not so large since the L- 
H-A contribution dominates p. For the 4 detector case, 
L-H-J-A network gives the largest value of piyr = 20.5. 
L-H-V-A network also gives similar results. In case of the 
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5 detector network, piyr = 21.4 and T^^^"^ = 7.21 days. 

IV. SUMMARY 

In this article we address the question of observing 
a hotspot of stochastic GW using the cross-correlation 
statistic. The idea is to restrict the statistic to a single 
or few pixels in the sky and target possible point stochas- 
tic sources. A possible source which we pick is the Virgo 
cluster which could be a rich bed of rotating neutron stars 
containing an estimated number of 10^^. Out of these the 
rotating neutron stars emitting GW which fall into the 
bandwidth of the advanced detectors are primarily the 
millisecond neutron stars. We assume that the distribu- 
tion of such neutron stars follows a bimodal distribution 
similar to that of the radio millisecond pulsars observed 
in our galaxy. We then see that with advanced detec- 
tors the observation time required to accumulate SNR 
^ 3 is about an order of an year if the average elliptic- 
ity of neutron stars is £ ^ 10~^. Several baselines have 
been considered as well as multiple baselines correspond- 
ing to networks of detectors. In these calculations, the 
baselines that stand out are the two LIGO detectors and 
the LIGO Livingston and a LIGO like detector at the 
AIGO site in Australia. These baselines have the best 
sensitivity, because for these baselines, the detectors are 
almost co-aligned. In such cases, the observation time 
required to achieve SNR ^ 3 is about 20 days if we as- 
sume e — 10~^. The future proposed Einstein Telescope 
can easily detect the hotspot. In fact, Einstein Telescope 



would be the only detector which can observe the Virgo 
cluster hotspot if e ~ 10~^. In that case, only one Ein- 
stein Telescope will be sufficient to detect the Virgo clus- 
ter by cross-correlating with other detectors like LIGOs, 
Virgo and LCGT. 

Besides the Virgo cluster, there could be other candi- 
dates for hotspots such as the Andromeda galaxy or our 
own galactic centre. Although in these cases, the number 
of sources contributing to the GW background may be 
smaller than the Virgo cluster, their distances are much 
smaller, which makes up for the overall strength of the 
stochastic sources. 
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